Use of in vitro methods combined with in silico analysis to identify potential skin sensitizers in the Tox21 10K compound library.

Introduction: Skin sensitization, which leads to allergic contact dermatitis, is a key toxicological endpoint with high occupational and consumer prevalence. This study optimized several in vitro assays listed in OECD skin sensitization test guidelines for use on a quantitative high-throughput screening (qHTS) platform and performed in silico model predictions to assess the skin sensitization potential of prioritized compounds from the Tox21 10K compound library. Methods: First, we screened the entire Tox21 10K compound library using a qHTS KeratinoSensTM (KS) assay and built a quantitative structure-activity relationship (QSAR) model based on the KS results. From the qHTS KS screening results, we prioritized 288 compounds to cover a wide range of structural chemotypes and tested them in the solid phase extraction-tandem mass spectrometry (SPE-MS/MS) direct peptide reactivity assay (DPRA), IL-8 homogeneous time-resolved fluorescence (HTRF) assay, CD86 and CD54 surface expression in THP1 cells, and predicted in silico sensitization potential using the OECD QSAR Toolbox (v4.5). Results: Interpreting tiered qHTS datasets using a defined approach showed the effectiveness and efficiency of in vitro methods. We selected structural chemotypes to present this diverse chemical collection and to explore previously unidentified structural contributions to sensitization potential. Discussion: Here, we provide a skin sensitization dataset of unprecedented size, along with associated tools, and analysis designed to support chemical assessments.

Fingerprinting Interactions between Proteins and Ligands for Facilitating Machine Learning in Drug Discovery.

Molecular recognition is fundamental in biology, underpinning intricate processes through specific protein-ligand interactions. This understanding is pivotal in drug discovery, yet traditional experimental methods face limitations in exploring the vast chemical space. Computational approaches, notably quantitative structure-activity/property relationship analysis, have gained prominence. Molecular fingerprints encode molecular structures and serve as property profiles, which are essential in drug discovery. While two-dimensional (2D) fingerprints are commonly used, three-dimensional (3D) structural interaction fingerprints offer enhanced structural features specific to target proteins. Machine learning models trained on interaction fingerprints enable precise binding prediction. Recent focus has shifted to structure-based predictive modeling, with machine-learning scoring functions excelling due to feature engineering guided by key interactions. Notably, 3D interaction fingerprints are gaining ground due to their robustness. Various structural interaction fingerprints have been developed and used in drug discovery, each with unique capabilities. This review recapitulates the developed structural interaction fingerprints and provides two case studies to illustrate the power of interaction fingerprint-driven machine learning. The first elucidates structure-activity relationships in ß2 adrenoceptor ligands, demonstrating the ability to differentiate agonists and antagonists. The second employs a retrosynthesis-based pre-trained molecular representation to predict protein-ligand dissociation rates, offering insights into binding kinetics. Despite remarkable progress, challenges persist in interpreting complex machine learning models built on 3D fingerprints, emphasizing the need for strategies to make predictions interpretable. Binding site plasticity and induced fit effects pose additional complexities. Interaction fingerprints are promising but require continued research to harness their full potential.

3D-Suspension culture platform for high throughput screening of neurotoxic chemicals using LUHMES dopaminergic neurons.

Three-dimensional (3D) cell culture in vitro promises to improve representation of neuron physiology in vivo. This inspired development of a 3D culture platform for LUHMES (Lund Human Mesencephalic) dopaminergic neurons for high-throughput screening (HTS) of chemicals for neurotoxicity. Three culture platforms, adhesion (2D-monolayer), 3D-suspension, and 3D-shaken, were compared to monitor mRNA expression of seven neuronal marker genes, DCX, DRD2, ENO2, NEUROD4, SYN1, TH, and TUBB3. These seven marker genes reached similar maxima in all three formats, with the two 3D platforms showing similar kinetics, whereas several markers peaked earlier in 2D adhesion compared to both 3D culture platforms. The differentiated LUHMES (dLUHMES) neurons treated with ziram, methylmercury or thiram dynamically increased expression of metallothionein biomarker genes MT1G, MT1E and MT2A at 6 h. These gene expression increases were generally more dynamic in 2D adhesion cultures than in 3D cultures, but were generally comparable between 3D-suspension and 3D-u plate (low binding) platforms. Finally, we adapted 3D-suspension culture of dLUHMES and neural stem cells to 1536 well plates with a HTS cytotoxicity assay. This HTS assay revealed that cytotoxicity IC50 values were not significantly different between adhesion and 3D-suspension platforms for 31 of 34 (91%) neurotoxicants tested, whereas IC50 values were significantly different for at least two toxicants. In summary, the 3D-suspension culture platform for LUHMES dopaminergic neurons supported full differentiation and reproducible assay results, enabling quantitative HTS (qHTS) for cytotoxicity in 1536 well format with a Robust Z' score of 0.68.

High-throughput screening of novel TFEB agonists in protecting against acetaminophen-induced liver injury in mice.

Macroautophagy (referred to as autophagy hereafter) is a major intracellular lysosomal degradation pathway that is responsible for the degradation of misfolded/damaged proteins and organelles. Previous studies showed that autophagy protects against acetaminophen (APAP)-induced injury (AILI) via selective removal of damaged mitochondria and APAP protein adducts. The lysosome is a critical organelle sitting at the end stage of autophagy for autophagic degradation via fusion with autophagosomes. In the present study, we showed that transcription factor EB (TFEB), a master transcription factor for lysosomal biogenesis, was impaired by APAP resulting in decreased lysosomal biogenesis in mouse livers. Genetic loss-of and gain-of function of hepatic TFEB exacerbated or protected against AILI, respectively. Mechanistically, overexpression of TFEB increased clearance of APAP protein adducts and mitochondria biogenesis as well as SQSTM1/p62-dependent non-canonical nuclear factor erythroid 2-related factor 2 (NRF2) activation to protect against AILI. We also performed an unbiased cell-based imaging high-throughput chemical screening on TFEB and identified a group of TFEB agonists. Among these agonists, salinomycin, an anticoccidial and antibacterial agent, activated TFEB and protected against AILI in mice. In conclusion, genetic and pharmacological activating TFEB may be a promising approach for protecting against AILI.

Exploring NCATS in-house biomedical data for evidence-based drug repurposing.

Drug repurposing is a strategy for identifying new uses of approved or investigational drugs that are outside the scope of the original medical indication. Even though many repurposed drugs have been found serendipitously in the past, the increasing availability of large volumes of biomedical data has enabled more systemic, data-driven approaches for drug candidate identification. At National Center of Advancing Translational Sciences (NCATS), we invent new methods to generate new data and information publicly available to spur innovation and scientific discovery. In this study, we aimed to explore and demonstrate biomedical data generated and collected via two NCATS research programs, the Toxicology in the 21st Century program (Tox21) and the Biomedical Data Translator (Translator) for the application of drug repurposing. These two programs provide complementary types of biomedical data from uncovering underlying biological mechanisms with bioassay screening data from Tox21 for chemical clustering, to enrich clustered chemicals with scientific evidence mined from the Translator towards drug repurposing. 129 chemical clusters have been generated and three of them have been further investigated for drug repurposing candidate identification, which is detailed as case studies.

Tautomeric cinnamoylphloroglucinol-monoterpene adducts from Cleistocalyx operculatus and their antiviral activities.

Guided by 1H NMR spectroscopic experiments using the characteristic enol proton signals as probes, three pairs of new tautomeric cinnamoylphloroglucinol-monoterpene adducts (1-3) were isolated from the buds of Cleistocalyx operculatus. Their structures with absolute configurations were established by spectroscopic analysis, modified Mosher's method, and quantum chemical electronic circular dichroism calculation. Compounds 1-3 represent a novel class of cinnamoylphloroglucinol-monoterpene adducts featuring an unusual C-4-C-1' linkage between 2,2,4-trimethyl-cinnamyl-ß-triketone and modified linear monoterpenoid motifs. Notably, compounds 1-3 exhibited significant in vitro antiviral activity against respiratory syncytial virus (RSV).

High-Throughput Screening to Advance In Vitro Toxicology: Accomplishments, Challenges, and Future Directions.

Traditionally, chemical toxicity is determined by in vivo animal studies, which are low throughput, expensive, and sometimes fail to predict compound toxicity in humans. Due to the increasing number of chemicals in use and the high rate of drug candidate failure due to toxicity, it is imperative to develop in vitro, high-throughput screening methods to determine toxicity. The Tox21 program, a unique research consortium of federal public health agencies, was established to address and identify toxicity concerns in a high-throughput, concentration-responsive manner using a battery of in vitro assays. In this article, we review the advancements in high-throughput robotic screening methodology and informatics processes to enable the generation of toxicological data, and their impact on the field; further, we discuss the future of assessing environmental toxicity utilizing efficient and scalable methods that better represent the corresponding biological and toxicodynamic processes in humans.

Auranofin targets UBA1 and enhances UBA1 activity by facilitating ubiquitin trans-thioesterification to E2 ubiquitin-conjugating enzymes.

UBA1 is the primary E1 ubiquitin-activating enzyme responsible for generation of activated ubiquitin required for ubiquitination, a process that regulates stability and function of numerous proteins. Decreased or insufficient ubiquitination can cause or drive aging and many diseases. Therefore, a small-molecule enhancing UBA1 activity could have broad therapeutic potential. Here we report that auranofin, a drug approved for the treatment of rheumatoid arthritis, is a potent UBA1 activity enhancer. Auranofin binds to the UBA1's ubiquitin fold domain and conjugates to Cys1039 residue. The binding enhances UBA1 interactions with at least 20 different E2 ubiquitin-conjugating enzymes, facilitating ubiquitin charging to E2 and increasing the activities of seven representative E3s in vitro. Auranofin promotes ubiquitination and degradation of misfolded ER proteins during ER-associated degradation in cells at low nanomolar concentrations. It also facilitates outer mitochondrial membrane-associated degradation. These findings suggest that auranofin can serve as a much-needed tool for UBA1 research and therapeutic exploration.

Exploring NCATS In-House Biomedical Data for Evidence-based Drug Repurposing.

Drug repurposing is a strategy for identifying new uses of approved or investigational drugs that are outside the scope of the original medical indication. Even though many repurposed drugs have been found serendipitously in the past, the increasing availability of large volumes of biomedical data has enabled more systemic, data-driven approaches for drug candidate identification. At National Center of Advancing Translational Sciences (NCATS), we invent new methods to generate new data and information publicly available to spur innovation and scientific discovery. In this study, we aimed to explore and demonstrate biomedical data generated and collected via two NCATS research programs, the Toxicology in the 21st Century program (Tox21) and the Biomedical Data Translator (Translator) for the application of drug repurposing. These two programs provide complementary types of biomedical data from uncovering underlying biological mechanisms with bioassay screening data from Tox21 for chemical clustering, to enrich clustered chemicals with scientific evidence mined from the Translator towards drug repurposing. 129 chemical clusters have been generated and three of them have been further investigated for drug repurposing candidate identification, which is detailed as case studies.

In silico modeling-based new alternative methods to predict drug and herb-induced liver injury: A review.

New approach methods (NAMs) have been developed to predict a wide range of toxicities through innovative technologies. Liver injury is one of the most extensively studied endpoints due to its severity and frequency, occurring among populations that consume drugs or dietary supplements. In this review, we focus on recent developments of in silico modeling for liver injury prediction using deep learning and in vitro data based on adverse outcome pathways (AOPs). Despite these models being mainly developed using datasets generated from drug-like molecules, they were also applied to the prediction of hepatotoxicity caused by herbal products. As deep learning has achieved great success in many different fields, advanced machine learning algorithms have been actively applied to improve the accuracy of in silico models. Additionally, the development of liver AOPs, combined with big data in toxicology, has been valuable in developing in silico models with enhanced predictive performance and interpretability. Specifically, one approach involves developing structure-based models for predicting molecular initiating events of liver AOPs, while others use in vitro data with structure information as model inputs for making predictions. Even though liver injury remains a difficult endpoint to predict, advancements in machine learning algorithms and the expansion of in vitro databases with relevant biological knowledge have made a huge impact on improving in silico modeling for drug-induced liver injury prediction.

In vitro profiling of pesticides within the Tox21 10K compound library for bioactivity and potential toxicity.

Pesticides include a diverse class of toxic chemicals, often having numerous modes of actions when used in agriculture against targeted organisms to control insect infestation, halt unwanted vegetation, and prevent the spread of disease. In this study, the in vitro assay activity of pesticides within the Tox21 10K compound library were examined. The assays in which pesticides showed significantly more activities than non-pesticide chemicals revealed potential targets and mechanisms of action for pesticides. Furthermore, pesticides that showed promiscuous activity against many targets and cytotoxicity were identified, which warrant further toxicological evaluation. Several pesticides were shown to require metabolic activation, demonstrating the importance of introducing metabolic capacity to in vitro assays. Overall, the activity profiles of pesticides highlighted in this study can contribute to the knowledge gaps surrounding pesticide mechanisms and to the better understanding of the on- and off-target organismal effects of pesticides.

Evaluating the utility of a high throughput thiol-containing fluorescent probe to screen for reactivity: A case study with the Tox21 library.

High-throughput screening (HTS) assays for bioactivity in the Tox21 program aim to evaluate an array of different biological targets and pathways, but a significant barrier to interpretation of these data is the lack of high-throughput screening (HTS) assays intended to identify non-specific reactive chemicals. This is an important aspect for prioritising chemicals to test in specific assays, identifying promiscuous chemicals based on their reactivity, as well as addressing hazards such as skin sensitisation which are not necessarily initiated by a receptor-mediated effect but act through a non-specific mechanism. Herein, a fluorescence-based HTS assay that allows the identification of thiol-reactive compounds was used to screen 7,872 unique chemicals in the Tox21 10K chemical library. Active chemicals were compared with profiling outcomes using structural alerts encoding electrophilic information. Random Forest classification models based on chemical fingerprints were developed to predict assay outcomes and evaluated through 10-fold stratified cross validation (CV). The mean CV Balanced Accuracy of the validation set was 0.648. The model developed shows promise as a tool to screen untested chemicals for their potential electrophilic reactivity based solely on chemical structural features.

Quantitative Bioactivity Signatures of Dietary Supplements and Natural Products.

Dietary supplements and natural products are often marketed as safe and effective alternatives to conventional drugs, but their safety and efficacy are not well regulated. To address the lack of scientific data in these areas, we assembled a collection of Dietary Supplements and Natural Products (DSNP), as well as Traditional Chinese Medicinal (TCM) plant extracts. These collections were then profiled in a series of in vitro high-throughput screening assays, including a liver cytochrome p450 enzyme panel, CAR/PXR signaling pathways, and P-glycoprotein (P-gp) transporter assay activities. This pipeline facilitated the interrogation of natural product-drug interaction (NaPDI) through prominent metabolizing pathways. In addition, we compared the activity profiles of the DSNP/TCM substances with those of an approved drug collection (the NCATS Pharmaceutical Collection or NPC). Many of the approved drugs have well-annotated mechanisms of action (MOAs), while the MOAs for most of the DSNP and TCM samples remain unknown. Based on the premise that compounds with similar activity profiles tend to share similar targets or MOA, we clustered the library activity profiles to identify overlap with the NPC to predict the MOAs of the DSNP/TCM substances. Our results suggest that many of these substances may have significant bioactivity and potential toxicity, and they provide a starting point for further research on their clinical relevance.

Three-Dimensional Structural Insights Have Revealed the Distinct Binding Interactions of Agonists, Partial Agonists, and Antagonists with the µ Opioid Receptor.

The United States is experiencing the most profound and devastating opioid crisis in history, with the number of deaths involving opioids, including prescription and illegal opioids, continuing to climb over the past two decades. This severe public health issue is difficult to combat as opioids remain a crucial treatment for pain, and at the same time, they are also highly addictive. Opioids act on the opioid receptor, which in turn activates its downstream signaling pathway that eventually leads to an analgesic effect. Among the four types of opioid receptors, the µ subtype is primarily responsible for the analgesic cascade. This review describes available 3D structures of the µ opioid receptor in the protein data bank and provides structural insights for the binding of agonists and antagonists to the receptor. Comparative analysis on the atomic details of the binding site in these structures was conducted and distinct binding interactions for agonists, partial agonists, and antagonists were observed. The findings in this article deepen our understanding of the ligand binding activity and shed some light on the development of novel opioid analgesics which may improve the risk benefit balance of existing opioids.

Identification of Potent and Selective Acetylcholinesterase/Butyrylcholinesterase Inhibitors by Virtual Screening.

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) play important roles in human neurodegenerative disorders such as Alzheimer's disease. In this study, machine learning methods were applied to develop quantitative structure-activity relationship models for the prediction of novel AChE and BChE inhibitors based on data from quantitative high-throughput screening assays. The models were used to virtually screen an in-house collection of ∼360K compounds. The optimal models achieved good performance with area under the receiver operating characteristic curve values ranging from 0.83 ± 0.03 to 0.87 ± 0.01 for the prediction of AChE/BChE inhibition activity and selectivity. Experimental validation showed that the best-performing models increased the assay hit rate by several folds. We identified 88 novel AChE and 126 novel BChE inhibitors, 25% (AChE) and 53% (BChE) of which showed potent inhibitory effects (IC50 < 5 µM). In addition, structure-activity relationship analysis of the BChE inhibitors revealed scaffolds for chemistry design and optimization. In conclusion, machine learning models were shown to efficiently identify potent and selective inhibitors against AChE and BChE and novel structural series for further design and development of potential therapeutics against neurodegenerative disorders.

In vivo quantitative high-throughput screening for drug discovery and comparative toxicology.

Quantitative high-throughput screening (qHTS) pharmacologically evaluates chemical libraries for therapeutic uses, toxicological risk and, increasingly, for academic probe discovery. Phenotypic high-throughput screening assays interrogate molecular pathways, often relying on cell culture systems, historically less focused on multicellular organisms. Caenorhabditis elegans has served as a eukaryotic model organism for human biology by virtue of genetic conservation and experimental tractability. Here, a paradigm enabling C. elegans qHTS using 384-well microtiter plate laser-scanning cytometry is described, in which GFP-expressing organisms revealing phenotype-modifying structure-activity relationships guide subsequent life-stage and proteomic analyses, and Escherichia coli bacterial ghosts, a non-replicating nutrient source, allow compound exposures over two life cycles, mitigating bacterial overgrowth complications. We demonstrate the method with libraries of anti-infective agents, or substances of toxicological concern. Each was tested in seven-point titration to assess the feasibility of nematode-based in vivo qHTS, and examples of follow-up strategies were provided to study organism-based chemotype selectivity and subsequent network perturbations with a physiological impact. We anticipate that this qHTS approach will enable analysis of C. elegans orthologous phenotypes of human pathologies to facilitate drug library profiling for a range of therapeutic indications.

Lestaurtinib induces DNA damage that is related to estrogen receptor activation.

A number of chemicals in the environment pose a threat to human health. Recent studies indicate estradiol induces DNA damage through the activation of the estrogen receptor alpha (ERα). Given that many environmental chemical compounds act like hormones once they enter the human body, it is possible that they induce DNA damage in the same way as estradiol, which is of great concern to females with the BRCA1 mutation. In this study, we developed an antibody-based high content method measuring γH2AX, a biomarker for DNA damage, to test a subset of 907 chemical compounds in MCF7 cells. The assay was optimized for a 1536 well plate format and had a satisfactory assay performance with Z-factor of 0.67. From the screening, we identified 128 compounds that induce γH2AX expression in the cells. These compounds were further examined for their γH2AX induction in the presence of an ER inhibitor, tamoxifen. After tamoxifen treatment, four compounds induced less γH2AX expression compared to those without tamoxifen treatment, suggesting these compounds induced γH2AX that is related to ERα activation. These four compounds were chosen for further studies to assess their ERα activating capability and c-MYC induction. Only lestaurtinib, a selective tyrosine kinase inhibitor, induced ERα activation, which was confirmed by both ERα beta-lactamase reporter gene assay and molecular docking analysis. Lestaurtinib also increased c-MYC expression, a target gene of ERα signaling, measured by the quantitative PCR method. This data suggests that lestaurtinib acts as a DNA damage inducer that is related to ERα activation.

Predictive Models for Human Cytochrome P450 3A7 Selective Inhibitors and Substrates.

Inappropriate use of prescription drugs is potentially more harmful in fetuses/neonates than in adults. Cytochrome P450 (CYP) 3A subfamily undergoes developmental changes in expression, such as a transition from CYP3A7 to CYP3A4 shortly after birth, which provides a potential way to distinguish medication effects on fetuses/neonates and adults. The purpose of this study was to build first-in-class predictive models for both inhibitors and substrates of CYP3A7/CYP3A4 using chemical structure analysis. Three metrics were used to evaluate model performance: area under the receiver operating characteristic curve (AUC-ROC), balanced accuracy (BA), and Matthews correlation coefficient (MCC). The performance varied for each CYP3A7/CYP3A4 inhibitor/substrate model depending on the data set type, model type, rebalancing method, and specific feature set. For the active inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.77 ± 0.01 to 0.84 ± 0.01. For the selective inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.72 ± 0.02 to 0.79 ± 0.04. The predictive power of the optimal models was validated by compounds with known potencies as CYP3A7/CYP3A4 inhibitors or substrates. In addition, we identified structural features significant for CYP3A7/CYP3A4 selective or common inhibitors and substrates. In summary, the top performing models can be further applied as a tool to rapidly evaluate the safety and efficacy of new drugs separately for fetuses/neonates and adults. The significant structural features could guide the design of new therapeutic drugs as well as aid in the optimization of existing medicine for fetuses/neonates.

Prediction of drug-induced liver injury and cardiotoxicity using chemical structure and in vitro assay data.

Drug-induced liver injury (DILI) and cardiotoxicity (DICT) are major adverse effects triggered by many clinically important drugs. To provide an alternative to in vivo toxicity testing, the U.S. Tox21 consortium has screened a collection of ∼10K compounds, including drugs in clinical use, against >70 cell-based assays in a quantitative high-throughput screening (qHTS) format. In this study, we compiled reference compound lists for DILI and DICT and compared the potential of Tox21 assay data with chemical structure information in building prediction models for human in vivo hepatotoxicity and cardiotoxicity. Models were built with four different machine learning algorithms (e.g., Random Forest, Naïve Bayes, eXtreme Gradient Boosting, and Support Vector Machine) and model performance was evaluated by calculating the area under the receiver operating characteristic curve (AUC-ROC). Chemical structure-based models showed reasonable predictive power for DILI (best AUC-ROC = 0.75 ± 0.03) and DICT (best AUC-ROC = 0.83 ± 0.03), while Tox21 assay data alone only showed better than random performance. DILI and DICT prediction models built using a combination of assay data and chemical structure information did not have a positive impact on model performance. The suboptimal predictive performance of the assay data is likely due to insufficient coverage of an adequately predictive number of toxicity mechanisms. The Tox21 consortium is currently expanding coverage of biological response space with additional assays that probe toxicologically important targets and under-represented pathways that may improve the prediction of in vivo toxicity such as DILI and DICT.

Use of Tox21 Screening Data to Evaluate the COVID-19 Drug Candidates for Their Potential Toxic Effects and Related Pathways.

Currently, various potential therapeutic agents for coronavirus disease-2019 (COVID-19), a global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), are being investigated worldwide mainly through the drug repurposing approach. Several anti-viral, anti-bacterial, anti-malarial, and anti-inflammatory drugs were employed in randomized trials and observational studies for developing new therapeutics for COVID-19. Although an increasing number of repurposed drugs have shown anti-SARS-CoV-2 activities in vitro, so far only remdesivir has been approved by the US FDA to treat COVID-19, and several other drugs approved for Emergency Use Authorization, including sotrovimab, tocilizumab, baricitinib, paxlovid, molnupiravir, and other potential strategies to develop safe and effective therapeutics for SARS-CoV-2 infection are still underway. Many drugs employed as anti-viral may exert unwanted side effects (i.e., toxicity) via unknown mechanisms. To quickly assess these drugs for their potential toxicological effects and mechanisms, we used the Tox21 in vitro assay datasets generated from screening ∼10,000 compounds consisting of approved drugs and environmental chemicals against multiple cellular targets and pathways. Here we summarize the toxicological profiles of small molecule drugs that are currently under clinical trials for the treatment of COVID-19 based on their in vitro activities against various targets and cellular signaling pathways.